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TECHNICAL NOTE NO . I 729 


FLIGHT DETERMINATION OF VING AND TAIL LOADS ON A FIGHTER - 

TYPE AIRPLANE BY MEANS OF STRAIN -GAGE MEASUREMENTS 
By VI Hi am S. Aiken, Jr. 


SUMMARY 


A flight investigation was conducted, to determine the contri- 
hutions of wing, tall, and fuselage to the total airplane lift of 
a propeller-driven fighter-type airplane . The tests covered a 
Mach number range from 0.2 to 0.8. The loads on the various airplane 
components were measured "by the use of calibrated strain-gage 
installations located at the roots of the wings and horizontal tail 
surfaces . 

The results of the investigation showed that the center of pres- 
sure of the wing additional air load remained constant near the 
lifting-line -theory value up to the highest Mach number (0.8) attained 
in the tests. The fuselage load -per unit normal acceleration appeared 
to he unaffected by Mach number changes. Adequate correlation was 
obtained between various tail-load parameters derived from the flight- 
test results and from wind-tunnel tests. Strain-gage measurements 
have been found to offer an accurate and convenient method for 
determining loads in flight. 


INTRODUCTION 


Although there are numerous wind-tunnel tests of airfoils and 
airpl ane models available at speeds which extend into the transonic 
range, the flight-loads data available for correlation are not only 
meager but for the most part are confined to relatively low Mach 
numbers. In order to obtain flight-loads data at high subsonic speeds, 
the National Advisory Committee for Aeronautics investigated a 
specially reinforced fighter-type airplane for the purpose of 
determining buffeting and transient loads . Some results of the 
buffeting tests have been published in reference 1. The strain- 
gage instrumentation used in the buffeting tests was adequate for 
the evaluation of wing additional-air-load center of pressure. 
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fuselage load per unit normal acceleration, wing-fuselage zero-lift 
pitching-moment coefficient, wing- fuselage aerodynamic center, and 
tail loads due to angular acceleration in pitch "below the buffeting 
boundary . The purposes of this paper are to present the results 
of the flight investigations "below the buffeting bo undar y; and to 
indicate the correlation existing between the quantities mentioned 
and the wind-tunnel-test results or theoretical calculations, or 
both- 


SYMBOLS 


B 


aerodynamic bending moment, inch-pounds 


Sh 


aerodynamic shear, pounds 


F 


fuselage load, pounds 


y C p spanwise distance of aerodynamic center of pressure 

Cram reference station, inches 

M Mach number 


1 

¥ 

X Y 

s 


g 


c 


N 




dynamic pressure, pounds per square foot 
airplane weight, pounds 

airplane moment of inertia in pitch, slug-feet square 
radius of gyration in pitch, feet 
wing area, square feet 

airplane normal acceleration at center of gravity 
(measured perpendicular to thrust line), g units 

acceleration due to gravity, feet per second per second 

airplane normal-force coefficient 



horizontal-tail load, pounds 

portion of horizontal-tail load associated with pitching 
acceleration, pounds 

pitching-momBnt coefficient of wing-fuselage combination 
at zero lift and M = 0 
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Cm o 

c 

X 

h 

9 angular acceleration in pitch, radians per second 

per second 

y lateral distance from airplane center line, feet 

t wing span, feet 

Subscripts : 

a additional air load 

-q. function of dynamic pressure 

L left wing 

R right wing 


wing section pitching-moment coefficient at zero lift 
mean aerodynamic chord, feet 

distance from wing-fuselage aerodynamic center to 
horizontal-tail center of pressure, positive 
forward, feet 

distance from airplane center of gravity to wing- 

fuselage aerodynamic center, positive forward, feet 


APPARATUS AND TESTS 
Airplane 


The airplane used in theBe tests was of conventional external 
appearance j however, the horizontal tail, fuselage, wing, canopy, 
and cowling had been heavily reinforced to provide an extra safety 
margin against structural failure in the investigation of buffeting 
loads. Pertinent geometric characteristics of the airplane are 
gLvefi in a three-view diagram, figure 1, and in table I. 


Instrumentation 

Standard NACA recording instruments were used to record time 
histories of impact -pres sure ; pressure altitude; normal, longitudinal, 
and lateral center-of- gravity accelerations; right and left aileron 
and elevator positions ; rate of pitch ; and elevator and aileron 
stick forces . A mechanical timer was used to correlate data from all 
recording instruments. 
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The airspeed head was mounted on a hoam extending 1.2 local 
chords ahead of the leading edge of the right wing tip . The 
NACA airspeed-altitude recorder was located in the right wing to 
minimize any lag effects • The entire airspeed system was calibrated 
for position error up to a Mach number of 0.78. 


Normal accelerations at the center line of the horizontal tail 
were measured "by means of a strain-gage-type accelerometer and re- 
corded hy a multiple recording oscillograph. 


The strains at various parts of the wing and horizontal tall 
structure were measured with wire -re si stance strain gages connected 
as four-arm. bridges and recorded hy a multiple recording oscillograph. 
Strain- gage bridges on the shear webs of wing and stabilizer spars were 
used to measure shear strains; strain-gage bridges on the spar flanges 
were used to measure bending strains . The strain-gage bridges were 


mounted at the 30-inch = ^*-*-35^ and 35~inch ~ 0*158) semispan 

stations on the right and left wings, respectively, and at the 
9 -inch semispan stations on the left and right stabilizers. In order 
to obtain calibrations for the conversion of measured strains into 
bending moments and shears, tests were carried out in the Langley air- 
craft loads calibration laboratory. The scope of the calibrations 
was such that the following accuracies for loads measurements could 
be estimated: 


Wing shears, pounds per side 1200 

Wing bending moments, inch -pounds per side ± 10,000 

Total horizontal- tail shear, pounds ±50 


Tests 

The tests reported herein all apply to load-factor conditions 
within the V-n envelope where no buffeting occurred. The data presented 
were obtained during flight at altitudes ranging from 5 >000 feet 
to 30,000 feet. In these tests the wing shears and bending moments 
were measured in either abrupt pull-ups to just below the stall or 
in a piloting sequence consisting of a rapid push-down followed by 
a rapid pull-up. An attempt was made to hold the Mach number constant 
and the ailerons neutral during these maneuvers. The average variation 

of Mach number in a given run was about ±jj? percent . 

Maneuvering tail loads were determined in slow turns as well as 
in the push-down - pull-up maneuvers . 

All flight tests were made with the airplane in the clean con- 
dition, with normal rated power, and with the center of gravity located 
at 25-1 percent mean aerodynamic chord. 
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RESULTS AND DISCUSSION 
Spanwiee Center of Pressure of Wing Load 


On the "basis of experience the lift or load distribution over 
a wing is usu all y considered to consist of two main components . 

One component, the basic distribution, consists of the distribution 
with a net lift of zero and is caused by geometric and elastic 
twist ^due to c^^, aileron deflection, an d rate of roll. This 

basic distribution is independent of angle of attack or normal-force 
coefficient. The other component, termed the "additional distribution," 
depends only on the angle of attack or normal-force coefficient at 
a given Mach number and dynamic pressure. 

With the load divided as stated, the bending moment at the wing 
root may be expressed by two terms as follows : 

B = + B a 

where B a is obtained from the equation 

B a = 


so that 


B B q. + ^cp^ 

In the foregoing equations the subscript q. refers to the bending- 
moment increment due to the basic air load, which varies with dynamic 
pressure, and the subscript a refers to the additional air loads. 

If during a given maneuver the dynamic pressure and the twist 
characteristics - that is, the twist, aileron deflection, and rate 
of roll - remain constant, then the center of pressure of the additional 
air load can be determined from the slope of the curve of bending 
moment plotted against shear. 

Fligdrb-test results .- in the flight bests, results from 3^- abrupt 
push-down and pull-up maneuvers were available in which the left- and 
right-wing loads were measured at pressure altitudes from 5,000 feet 
to 30,000 feet and at Mach numbers from 0.2 to 0.8. The maneuvers 
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covered, a noxmal-force-coefficient range from -0.66 to 1.25* A 
typical maneuver illustrating the type of data available for analysis 
is shown in figure 2 for an abrupt push-down - pull-up maneuver 
executed at 10,000 feet. It will he noted in this figure that the 
Mach number is reasonably constant throughout the maneuver because 
the maneuver was abrupt . There was less than a 100-foot change in 
pressure altitude. 

The strain gages measure a structural load, that is, an aero- 
dynamic load plus an inertia load; therefore, acceleration data 
were used to correct the measured loads for inertia effects by 
simply adding values of the wing weight multiplied by n - 1 to 
the shear measurements outboard of the strain- gage measuring station. 

A similar correction was applied to be ndin g-moment results. Now, 
as stated previously, if the dynamic pressure remains constant during 
a maneuver, the center of pressure of the additional air load is 
determined by taking the slope of the curve of bending moment plotted 
against shear. The aerodynamic bending moments for the left and right 
wings are plotted against the corresponding aerodynamic shears in 
figure 3 for the maneuver illustrated in figure 2 . The slope of each 
line was found by a least- squares method and the centers of pressure 
thus determined were 75-9 inches 'for the left wing and 79*8 incheB 
for the right wing. These centers of pressure are for the load outboard 
of the strain-gage semispan station, but an approximate check indicates 
that the centers of pressure for these loads on each wing are roughly 
at the same spanwise location as follows : 

y 0 -n = 75*9 + 35 = 110-9 inches 

cp L 

y = 79.8 + 30 = 109*8 inches 
Cp R 


The procedure outlined was followed for all the maneuvers 
considered and the results, in terms of the center of pressure of the 
additional air load outboard of the strain-gage station, are given 
in figure 4 as a function of Mach number. The theoretical rigid-wing 
additional-air-load centers of pressure for the test airplane were 
computed by lifting- line and strip theories and are shown in figure k 
as the solid and dashed lines, respectively. The experimental data 
show no trend for these centers of pressure to move either outboard 
or inboard with Mach number but rather tend to scatter about the line 
predicted by lifting- line theory. Sixty percent of the points shown 
in figure 4 represent deviations from lifting- line -theory values of 
less than the probable error associated with the determination of the 
experimental center of pressure . The average probable error for all 
points shown is ±1.63 inches. 


Wind-tunnel results .- An analysis of data from pressure-distribution 
measurements made on a "scale model of a prototype of the test airplane 
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was made to determine the variation of addltlonal-air-load center of 
pressure -with Mach number. The wing plan forms and aerodynamic charac- 
teristics of the model and the test airplane were essentially the same . 

y 

Pressure distributions were measured at semispan station — ^ equal 

to 0.211, 0.511 , 0.673.> 0*839j aid O.935. The tests were conducted in 
the Ames l6-foot high-speed wind tunnel to find the effect of aileron 
deflection on span loading. The additional-air-load centers of pressure 
for the model were calculated by using points at various normal- force 
coefficients at a constant Mach number and then deter minin g the slopes 
of the curves of bending moment plotted against shear. The results are 

shown in figure 5 Tor aileron angles of 9° and “9° • It should be re- 
membered that the additional air load is independent of ailBron angle . 
The centers of pressure are expressed in terms of the dimensions for 
the full-size airplane . There is an apparent outboard movement of the 
center of pressure in the Mach number range from 0-70 to 0.80, but 
this trend is not evident from flight- test results (fig. 4) and is 
somewhat questionable when the scatter of the data at constant Mach 
number for the two different aileron angles is considered, since the 
additional- air-load center of pressure is independent of aileron 
angle. The discrepancy may be due to differences in elastic charac- 
teristics of the model wing and the airplane wing. 

Plight data from -pressure-distribution measurements on two 
airplanes .- Although flight data on wing span loading are generally 
meager, tests on two airplanes (references 2 and 3) furnished a 
limited amount of data from which the additional-air-load center 
of pressure could be computed. The variation of additional-air-load 
center of pressure with Mach number is shown in figure 6(a) for the 
airplane of reference 2 and in figure 6(b) for the airplane of 
reference 3* Por the airplane in figure 6(a) there is a trend 
toward inboard movement of the center of pressure with increasing 
Mach number, but considerable error is possible in the determination 
of the points at M = 0 .70 and M = O.78 since only two airplane 
normal-force coefficients were available for the computation of the 
center of pressure. The data of figure 6(b) show no change in 
center of pressure with Mach number . The results shown in figures 4 
to 6, in general, indicate that the additional-air-load center 
of pressure of an unstalled unswept wing may be expected to remain 
constant up to a Mach number of O.85 and to be in good agreement with 
the lifting-line-theory values. 

Pusolage Loads 


The determination of fuselage loads by using strain-gage measurements 
as has been done in the case of the test airplane is relatively simple. 
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If (the normal acceleration, the airplane weight, the horizontal- 
tall load, and the left- and the right-wing shears are known at al 1 
times, the fuselage load can he found from the agnation 



The fuselage load thus includes the smaT 1 portion of the wing area in- 
hoard of the strain-gage stations . In order to determine the variation 
of fuselage load with Mach number, the push-down - pull-up maneuvers 
described previously were used. The fuselage load E was plotted 

against the normal acceleration; the slope of this line ^ is the 

fuselage load per uni t normal acceleration for a particular Mach number. 
A typical variation of fuselage load with normal acceleration is 
illustrated in figure 7 for a Mach number of 0.51 at a pressure 
altitude of 20,000 feet. The fuselage load per unit normal acceleration 
found from curves of this type is plotted in figure 8 as a function of 
Mach number. A theoretical value (1713 lh/g), computed by assuming that 
the load carried by the fuselage was equivalent to the load on the wing 
area inboard of the strain-gage stations and through the fuselage, 
passes approximately through the center of the data. No apparent 
variation of fuselage load with Mach number may he noted. , 


Horizontal-Tail loads 

The aerodynamic tail load may he expressed as 


L, = -C 
t m 


'0 \Ji - M 2 *t *t 


qSc ^ eNx 

*t 


( 1 ) 


The factor l/vl - M^ is introduced in this equation because the 
Cjjj^-term cannot otherwise he determined. The Mach number range is 

l imi ted to 0.6, and in this range the Glauert factor should h 9 adequate. 
In order to supplement data obtained in high-speed dives a series of 
very slow turns were made in which 9 was approximately equal to zero . 
Equation (l) may then he written as 



= -C 


Sc n Wx 
“0 2 t + q 


\[l - M 2 
st 
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or 




i/yi- 


M 


_ = _o -So „ Sx 
2 “o ^ 


\ll - # 


The measured values of 






In the Mach number range from 


q./ VI “ M 

0.2 to 0.8 'were plotted against Cjj\/l - M 2 as shown in figure 9* 

The data of figure 9 are for an airplane weight of 8750 pounds and 
a cent er-of -gravity position of 25*1 percent mean aerodynamic chord. 

The relative lack of scatter for this data Justifies the inclusion 
of the Glauert factor. A curve faired through this data is linear 
up to a normal-force coefficient of about 0 .5 * In this linear region 
the aerodynamic-center position x may be determined from the slope of 
the cupve as 




L / \l 1 - *£/ _ Jk = 


( C N n/* 7 ^) 


£& _ Sx 


Z t + x 


With x determined, C-^ may be found from the intercept of 


the curve at C. 


»P = 0 


as 


%) 


Z t + * 
Sc 


■/n/ 1 - ^ 


The variation of the quantity 


J t- 


L /\/l - M 2 


with 




as 


det'ermined 


from tests of a --scale model of an airplane (reference with the 

same basic design as the airplane shown in figure 1, is shown 
in figure 9 • The agreement in the lower normal-f orce-coefficient range 
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is considered good. In. table H are listed the tail-load parameters 
which may he derived from the curves shown in figure 9 • The agreement • 
between the low-speed values of for the flight-test results and 

wind-tunnel results is excellent as is the agreement between the 
tail-off aerodynamic -center position x and the tail load per unit 

normal acceleration at low values of for both flight and 

wind-tunnel results. The values of tail load per unit normal acceleration 

listed in table H apply only at the listed values of “bn/ 1 - since 
actual tail loads would have to be determined by an integration process 

at values of C J 1 above 0 . 4 . 


The variation of the aerodynamic -center position in percent of the 
mean aerodynamic chord from the leading edge is given in figure 10 as a 


function of 


°H\/l 


for the test airplane . The values in this 


figure were obtained from the slope of the line through the experi- 
mental points in figure 9. The tail-off aerodynamic center shifted 
forward with increase in airplane normal-force coefficient. 


The determination of aerodynamic tail loads due to angular 
acceleration in pitch must be accomplished indirectly by subtracting 
calculated values of balancing tail load from measured tail loads 
and plotting the resulting differences against airplane angular 
acceleration. The results of these computations made for several 
abrupt push-down - pull-up maneuvers are shown in figure 11 . The 
angular accelerations were determined from normal- acceleration 
measurements at the center of gravity and the tail of the test air- 
plane. The theoretical values of were calculated from the 

expression 


VK 2 .. 
H- e - = — 9 


SH 


where ¥ equals 8j 50 pounds, ky equals 5-5 feet, g equals 32.2 feet 

per Lecond per second, and x, equals -16.78 feet. For these values 
dl t -- 

— equals -490 pounds per radian per second per second. The 

d0 

agreement of the experimental values with the calculated values is 
considered adequate since the determination of the experimental values 
required the use of smal 1 differences of large numbers. 
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CONCLUSIONS " 


From the results of the flight tests and correlation -with some 
wind- tunnel- test results, the following conclusions were indicated: 

1 . The results of the determination of wing additional- air-load 
center of pressure from strain-gage measurements of root tending 
mamenbs shears for the test airplane indicated good agreement 
with lif t ing - lin e theory and showed that the center of pressure of 
the additional air load might "be ejected to remain constant up to 

a Mach number of 0 . 8 . 

2. A comparison of wind-tunnel additional-air- load center-of- 
pressure results for a model of a prototype of the test airplane with 
flight results showed that the outboard shift of load predicted from 
wind-tunnel tests was not evident in flight. 

3- With the data determined from flight pressure-distribution 
measurements on the airplanes of NACA AKR No. hFD9 and NACA TN No. 11 8 l, 
it appeared that the additional-air-load center of pressure might rema i n 
constant to Mach numbers of O. 85 . 

The fuselage load per unit normal acceleration for the test 
airplane remained approximately constant in the Mach number range of 
the tests with the data tending to scatter about a theoretical value 
which was computed by assuming that the load carried by the fuselage 
was equivalent to the load on the wing area inboard of the strain-gage 
stations and through the fuselage . 

5 . The tail-load measurements for the test airplane showed 
excellent agreement between flight and wind-tunnel values of the 
zero-lift pitching-moment coefficient and between flight and wind- 
tunnel values of the tail-off aerodynamic-center position at low 
airplane normal-force coefficients <» 

6 . Flight tail-load measurements indicated that as the airplane 
normal-force coefficient increased there was a forward shift in 

the tail-off aerodynamic -center position. 

7 . Values of tail loads due to angular accelerations in pitch 
derived from flight measurements gave reasonable agreement with 
calculated values. 
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8. Because of the agreement Between the flight measurements and 
wind-tunnel or theoretical results, the use of properly calibrated 
strain-gage installations to measure wing and tail loads offered an 
accurate and convenient means of dete rmining over-all loads in flight . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., August 13, 19^8 
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TABLE I 

CHARACTEBISTICS OF AIRPLAHE 

Wing: 

Span, ft 37 *03 

Area, sq. ft 2^0.1 

Mean aerodynamlo chord, ft 6.63 

Airfoil Low drag 

Aspect ratio 5*7 

Horizontal tail: ■ 

Span, ft 13.18 

Area, sq. ft .4-1.0 

Incidence, cleg 1 

Weight during tests, lb 8750 

Center-of -gravity position during tests, 

percent M .A .C . * • 25 .1 

Approximate radius of gyration in pitch, ft ' 5*5 

Tail length, Z.J. , ft -15.84- 



TABLE! H 


COMPARISON OF FLIGHT AND WIND-TUNNEL TAIL-LOAD PARAMETERS 
[Center of gravity at 25 .1 percent M.A.C.J W « &750 -lb] 


Tail-load parameter 

Flight 

results 

Wind-tunnel 
results frcm 
reference 4 

Zero-lift tail-off pitoMng-mcanent coefficient, 0^ 

- 0.0520 

- 0.0528 

Tail-off aerodynamic -center position, distance from center of gravity 
to aerodynamio center, x, percent M.A.C..: 



At O'Jl - $ - 0 no 0.4 

- 10 .98 

- 10.60 

At 0-j.Wl - MP « 0.5 

-11,25 

- 10.60 

At C N ^1 - Mp = 0.7 

-13-17 

- 10.60 

At 0^1 - M 2 » 0.9 

- 25 .10 

- 10.60 

Tail load per unit normal acceleration, lb/g: 



At 0^1 - = 0 to 0.4 

380 

367 

At C^/l - # - 0,5 * - • • 

389 

367 

At C^l - - 0.7 

450 

367 

At Ojj^l - *£ » 0.9 

702 

367 








Three ~v/ew d/agram o-f Test a/rp/are. 





Perodj/ra/r/c shear, //> Perodyna/rlc bending /7xmeaf,ln.-tt> 
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^NACA. 



pfac/y rumfsr- 

f/gure & - t&r/of/o/y of fi/e»/cype /oac/pe/' c//?/f pormo/ a cce/erat/or? 
m/ff /ffaaf numter forte of o/spf?/?e. • 
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F/gure v9 —/ar/crf/on of fa//- /oad para/nefer wfh a/rp/aoe norma/ force 

coefficient for feof o/r/o/ane. d/go/one /re/gfi/ d?73& /aoi/nds- 
c enter- of- grar/fp /aos/f/or? ) /oercenf /neon aerodg'no/n/c 

c fiord. 






Angular acceleration, radians per sec z 

Figure // Variation of angular- acceleration fat! toad mth angular accelerator as detrr- 

m/ned from flight data for test airplane. 
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